Ab Initio Structural Energetics of Si3N4 Surfaces 
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Motivated by recent electron microscopy studies on the Si3N4/rare-earth oxide interfaces, the 
atomic and electronic structures of bare /3— Si3N4 surfaces are investigated from first principles. 
The equilibrium shape of a Si3N4 crystal is found to have a hexagonal cross section and a faceted 
dome-like base in agreement with experimental observations. The large atomic relaxations on the 
prismatic planes are driven by the tendency of Si to saturate its dangling bonds, which gives rise 
to resonant-bond configurations or planar sp^-type bonding. We predict three bare surfaces with 
lower energies than the open-ring (1010) surface observed at the interface, which indicate that 
non-stoichiometry and the presence of the rare-earth oxide play crucial roles in determining the 
termination of the Si3N4 matrix grains. 



PACS numbers: 68.35.Bs, 68.35.Md, 68.37.Lp 

The desirable mechanical and physical properties of 
silicon nitride ceramics Q in many high temperature ap- 
plications are hindered by their intrinsic brittleness which 
limits their wide-spread use and reliability as structural 
components. It has been empirically known for some 
time that this problem can be overcome by microstruc- 
tural and compositional design with sintering additives, 
in particular rare-earth oxides 2] . The resulting ceramic 
microstructure consists of elongated Si3N4 matrix-grains 
embedded in an intcrgranular, typically amorphous, rare- 
earth oxide phase. However, precise information about 
the structure and chemistry of the interface has been 
lacking for many years. Recently, there have been three 
experimental studies using scanning transmission elec- 
tron microscopy (STEM), which revealed for the first 
time important information about the atomic structures 
and bonding characteristics at the /3— Si3N4/rare-earth 
oxide interfaces 0. This exciting development provides 
a timely motivation for systematic theoretical studies of 
the interface, which should complement and aid in the 
interpretation of these experiments. 

Existing theoretical calculations in this field 4] , which 
have typically employed tight-binding, pair potential, 
and molecular orbital methods (a few of which are from 
first principles), have mainly focused on the bonding sites 
of additive atoms at the interface and the resulting elec- 
tronic structure. While these studies have provided im- 
portant information on the Si3N4/rare-earth oxide inter- 
face, our approach is to start with detailed first principles 
calculations on bare Si3N4 surfaces, which have not yet 
been performed. We believe that such an approach is 
the first natural step in a systematic understanding of 



the interface at the microscopic level. As it has been 
historically the case in ab initio investigations of several 
other important interfaces, the insights gained from this 
step can not only be used as a basis to elucidate the 
bonding characteristics of the additives, but also raise 
new questions which need to be addressed to have an 
in-depth understanding of the interface. In this Letter, 
we present, for the first time, results from first princi- 
ples calculations on low-index /3— Si3N4 surfaces focusing 
in particular on (i) the equilibrium shape of the crys- 
tal, and (ii) the atomic structure and stoichiometry at 
the prismatic plane (1010) surface, which is the relevant 
surface studied in recent STEM experiments [Fig. 1(a)]. 

Our calculations were performed within density func- 
tional theory using the projector augmented wave 
method 0. For exchange-correlation, we used the 
Perdew-Wang parametrization of the generalized gradi- 
ent approximation (GGA), and repeated some of the cal- 
culations using the Ceperley- Alder functional within the 
local density approximation (LDA) for comparison. An 
energy cutoff of 270 eV was used in all calculations. For 
structural optimization of the bulk, a F-centered 3x3x8 
k-point grid was employed. Doubling the k-point grid 
and increasing the cutoff to 400 eV had no apprecia- 
ble effect on the calculated structural parameters. The 
hexagonal unit cell of /3-Si3N4 contains 14 atoms with 
6 structural parameters (a, c, and 4 internal parame- 
ters). Half of the atoms are in the plane z = jc, and the 
other half are in the z — jc plane, as shown in Figure 
1(b). Si atoms are four-fold coordinated with N atoms 
in a slightly distorted tetrahedral configuration, and N 
atoms are threefold coordinated with Si atoms. The 
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FIG. 1: (a) Z-contrast image of a Si3N4 grain with a hexag- 
onal edge showing the interface with LU2O3. The grain is 
oriented with the [0001] projection (note the superimposed 
atomic structure), so that the (1010) prismatic boundary 
planes are set at an "edge-on" condition. Lu atoms are the 
bright spots attached in pairs at the termination of the hexag- 
onal open rings, (b) A perspective view of the /3— Si3N4 unit 
cell. Si and N atoms are shown in white and gray circles, 
respectively, (c) [0001] projected view of the unit cell. 



six structural parameters (a, c, ccsi, ysi, Xf^, Un) optimized 
with GGA (7.667 A, 2.928 A, 0.1752, 0.7692, 0.3299, 
0.031) and LDA (7.585 A, 2.895 A, 0.1738, 0.7675, 0.3301, 
0.0295) are in good agreement with previous calculations 
and experimental values of (7.608 A, 2.911 A, 
0.1733, 0.7694, 0.3323, 0.0314). The surface calculations 
were performed with the relaxed structural parameters 
using a slab geometry. Depending on the surface, we 
used hexagonal and simple or base-centered monoclinic 
supercells (Table I). For each surface, we considered up to 
4 different stoichiometric terminations to find the lowest- 
energy atomic configuration. We performed several tests 
to assess the convergence of the surface energy with re- 
spect to number of layers (up to 7) and the size of the 
vacuum region (up to 10 A). Although the convergence 
was observed to be rapid [with the exception of the (1120) 
surface], we used 5-layer slabs [10-layer for (1120)] with 
8 — 10 A vacuum. 

Experimentally, it is observed that the Si3N4 mi- 
crostructure consists of elongated structures with hexag- 
onal cross sections. The preferential growth along the 
c— axis is rather rapid, while the growth of the (1010) 
and (1120) prismatic planes is reaction limited and very 
sensitive to the type of additives 0. In order to under- 
stand the growth process of bare Si3N4 grains, we calcu- 
lated the surface energies of five low-index Si3N4 surfaces 
(Table I). The equilibrium shape of /3— Si3N4 calculated 
from these energies by the Wulff construction ^ is shown 
in Figure 2. The exposed faces consist of the {1120}, 
{1011}, {1121}, and {0001} families of planes making up 
~ 53 %, 25 %, 19 %, and 3 % of the total surface area, re- 
spectively. The calculated aspect ratio for bare surfaces 
between the crystal length and its width is 1.4, which 
is in agreement with the measured aspect ratios near 2 
for Si3N4 samples sintered without or only the minimum 



TABLE I; Supercell shape, parameters, and surface energies 
for the five lowest-index /3— Si3N4 surfaces. In the second 
column, the numbers in parentheses are the surface unit cell 
dimensions; M (monoclinic) and BCM (base-centered mono- 
clinic) refer to the lattice type of the supercell, and /3 is the 
monoclinic angle. The energies for (1120), (1121), and (1010) 
refer to the lowest-energy stoichiometric terminations, while 
the energy for (1010) is that of the open-ring surface. 
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FIG. 2: Two views of the equilibrium shape of a /9-Si3N4 
crystal calculated using the Wulff construction. 4 families 
of surfaces are exposed: {1120}, {0001}, {1121} (the shaded 
family of surfaces), and {1011} (the unshaded family of sur- 
faces surrounding the (0001) surface). The bold numbers in 
parentheses are the calculated surface energies in J/m^ The 
calculated L/W aspect ratio is 1.4. 



amount of the additive. We expect that in the presence 
of the rare-earth oxide additive, the energies of (1120) 
or (1010) will be lower, thereby increasing the calculated 
aspect ratio. We also note that our calculations predict 
an equilibrium shape, in which the base is not flat, but 
like a faceted dome even when only two tilted families 
of planes, {1121} and {1011}, are included. This is in 
agreement with the experimental observation on macro- 
scopic Si3N4 matrix-grains embedded in the intergranular 
phase, where the basal planes look atomically rough ic|]. 

It is particularly interesting to notice that of the two 
lowest-index prismatic planes rotated by 30° with respect 
to each other, it is the (1120) surface which is exposed in 
the Wulff construction and not the (1010) surface with 
open hexagonal rings, which has been shown to exhibit 
an abrupt interface with the rare-earth oxide additive in 
recent STEM experiments. The reason for the signifi- 
cantly lower surface energy of (1120) (1.95 J/m^) com- 
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pared to (1010) (2.57 J/m^) can be understood from the 
nature of the atomic relaxations and the way the dan- 
gling bonds of Si are saturated. On the ideal (unrelaxed) 
(1010) surface, there are one Si and one N atom, which 
have one dangling bond each. These are shown by Si2 
and Nl in Fig 3(a). When the surface is relaxed, while 
Nl still remains with a dangling bond, Si2 undergoes a 
considerable displacement ('^ 0.8 A) forming a new Si-Si 
bond with Si4 at 2.58 A, only 10 % larger than the bulk 
Si-Si distance [Fig. 3(b)]. As such, the dangling bond 
of Si2 is saturated, and the 5-fold coordinated Si4 atom 
has two resonant bonds, reminiscent of overcoordination 
of Si in negatively charged center defect The 
surface is still somewhat rough due to the presence of 
the Sil-Nl unit, which does not relax significantly. On 
the ideal (1120) surface, on the other hand, there are 
2 Si and 2 N atoms with one dangling bond each, la- 
beled by Si6, Si4, N5, and N3 in Fig. 3(c). When this 
surface is relaxed, the resulting atomic displacements for 
most of the atoms are very large. For example, the dis- 
placements of N7, Si5, and Si6 are 1.51 A, 1.49 A, and 
1.33 A, respectively. Such large relaxations saturate all 
the dangling bonds except for N3, resulting in a rather 
smooth surface with 7-fold, 3-fold, and 4-fold Si-N rings 
[Fig. 3(d)]. The only undercoordinated Si atom on the 
relaxed surface (Si5 with no dangling bonds on the ideal 
surface) exhibits an sp^— type bonding in an almost pla- 
nar coordination with 3 N atoms. 

We also investigated the relative structural stability of 
the (1010) bare surface as a function of stoichiometry. 
The formation energy of a surface, -Eform, can be written 
as [12 

£^form = ^^sTab " "SiMSi " m fJ-N , (1) 

where i?*°ab is the total energy of the slab, Ui is the num- 
ber of atoms of type i in the slab, and /^i's are the corre- 
sponding chemical potentials, which satisfy 3/zsi + 4/iN — 
MSiaNi.buik- Figure 4 shows the calculated energies as a 
function of the stoichiometry of the surface. All ener- 
gies are relative to the stoichiometric surface with open 
hexagonal rings, which we will refer to as the "open-ring 
surface" . This reference configuration is quite stable in 
the full allowed region of /iN with the exception of two 
other terminations. The first one, which we will refer to 
as the "half-surface" , is obtained by removing half the 
atoms from the open-ring surface as shown by the dot- 
ted line in Fig. 3(a). This surface, while stoichiometri- 
cally terminated and quite smooth even in its ideal struc- 
ture, does not have the open hexagonal rings, which are 
present at the Si3N4 rare-earth oxide interfaces studied 
in recent STEM experiments. When relaxed, the half- 
surface becomes 0.17 eV (GGA) and 0.1 eV (LDA) lower 
in energy compared to the open-ring surface. The only 
significant relaxation occurs for the Si atom (Si4) with 
the dangling bond which moves inward by ~ 0.64 A. As 
a result, the Si atom exhibits a sp^— type bonding in a 
nearly planar coordination with 3 N atoms [Fig. 3(e)]. 
One possibility for the discrepancy between the theoret- 




FIG. 3: [0001] projected, double-unit-cell views of (a) 2-layer 
ideal (1010) open-ring surface, (b) the top-layer of the re- 
laxed open-ring surface, (c) 2-layer ideal (1120) surface, (d) 
2-layer relaxed (1120) surface, (e) top-layer of the relaxed 
half-surface, and (f) relaxed structure of a non-stoichiometric 
(1010) surface with a missing SiN2 unit (1.5 layers shown). 
In (a), the index 1, for example, inside a white circle denotes 
the Sil atom mentioned in the text. If the atoms above the 
dashed line are removed, the ideal half-surface is obtained. 
If the atoms inside the dashed ellipse (Nl, Sil, and N4) are 
removed, the ideal structure of the non-stoichiometric surface 
(mentioned in the text) is obtained. In (c), the dashed region 
corresponds to a single unit cell of the (1120) surface. 



ical prediction of the half-surface having a lower energy 
than the experimentally observed open-ring surface at 
the rare-earth oxide interface is that theory incorrectly 
predicts the lowest energy termination of the bare (1010) 
Si3N4 surface. This could be verified by studying the 
bare surface with low-energy electron diffraction experi- 
ments. A more likely reason for the discrepancy, which 
would also explain the prediction of the (1120) surface 
as the lower-energy prismatic plane contrary to experi- 
mental observations, is that the rare-earth oxide additive 
changes the relative stability of different Si3N4 surfaces 
and different terminations of the same-index surfaces. 

The second termination, which results in a lower en- 
ergy surface under N-poor conditions, is obtained by re- 
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FIG. 4: (Color) The formation energies of various termina- 
tions of the (1010) surface relative to the open-ring surface 
as a function of stoichiometry Afi — fi]^ — |^tN2- (+) and 
(— ) signs correspond to taking out and adding the particular 
atom, respectively, from the open-ring surface shown in Fig. 
3(a) with all the atom indices. 

moving a SiN2 unit (Nl, N4, Sil) from the open-ring sur- 
face, as shown by the dashed ellipse in Fig. 3(a). In its 
ideal structure, this surface has 2 Si atoms, Si4 and Si2, 
with 1 and 2 dangling bonds, respectively, along with a 
N atom (N3) with one dangling bond. In spite of a large 
number of broken bonds, when the surface is relaxed, 
large atomic displacements of Si2 (1.29 A) and Si4 (1.27 
A) saturate all Si dangling bonds (via the formation of 
a new Si-Si bond at 2.5 A), which reduces the surface 
energy considerably. The relaxed surface is smooth and 
has 7-fold and 3-fold rings, similar to those observed on 
the (1120) surface [Fig. 3(f)]. It is interesting to note 
that under extreme N-poor conditions, the energy of this 



surface is 1.85 J/m^, which is lower than that of (1120), 
providing a possible explanation for the dominant obser- 
vation of (1010) surfaces in the STEM experiments. 

In summary, we have presented results from ab initio 
calculations on the atomic and electronic structures of 
bare /?— Si3N4 surfaces with emphasis on the prismatic 
plane (1010) surface observed in recent STEM experi- 
ments. The equilibrium shape of a macroscopic Si3N4 
crystal is found to have a hexagonal cross section, a 
faceted dome-like base, and an aspect ratio of 1.4, in 
agreement with experimental observations. We find large 
distortions on the prismatic planes driven primarily by 
the tendency of Si atoms to saturate their dangling bonds 
and achieve either resonant-bond, sp^ — , or sp^— bonded 
configurations. The stoichiometric (1120) surface, the 
(1010) "half-surface", and a non-stoichiometric (1010) 
surface obtained by removing a SiN2 unit are predicted 
to have lower energies than the open-ring (1010) surface. 
In light of the consistent experimental observations of the 
open-ring (1010) surface at the interface, the present re- 
sults obtained with state-of-the-art ab initio techniques 
strongly indicate that (i) the rare-earth oxide additive 
changes the relative stability of /3— Si3N4 surfaces, and 
(ii) non-stoichiometry (especially resulting from N-poor 
conditions) should play an important role in determining 
the termination of the Si3N4 matrix grains. We expect 
that future ab initio studies of the interface will not only 
focus on the bonding nature of the additives to the open- 
ring (1010) surface, but also address why the proposed 
mechanisms do not favor promotion of the (1120) ter- 
mination, which has a lower bare surface energy. This 
work was supported by DOE under grants No. DE- 
FG02-03ER15488 (JCI and SO), DE-AC05-03ER46057 
(AZ and NDB), and DE-AC02-05CH11231 (ROR), and 
by the ACS Petroleum Research Fund under grant No. 
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